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The recent discovery of a ghost-free, non-linear extension of the Fierz-Pauli theory of 
massive gravity, and its bigravity formulation, introduced new possibilities of interpreting 
cosmological observations, in particular, the apparent late-time accelerated expansion of 
the Universe. Here we discuss such possibilities by studying the background cosmology 
of the model and comparing its predictions to different cosmological measurements. We 
place constraints on the model parameters through an extensive statistical analysis of 
the model, and compare its viability to that of the standard model of cosmology. We 
demonstrate that the model can yield perfect fits to the data and is capable of explaining 
the cosmic acceleration in the absence of an explicit cosmological constant or dark energy, 
but there are a few caveats that must be taken into account when interpreting the results. 
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Recent work^"'' has solved a long-standing and highly- challenging puzzle in theoret- 
ical physics, that was first proposed by Fierz and Pauli in 1939:^ "does a consistent 
non-linear theory of massive gravity exist?" In Refs. 1 and 2, de Rham, Gabadadze 
and Tolly showed for the first time that it is possible to give gravity a mass while 
the theory remains free of the so-called Boulware and Deser ghosts^ (for a review 
of massive gravity, see Ref. 10). Later, in Refs. 3, 4, 5, 6 and 7, Hassan and Rosen 
demonstrated that such a ghost-free theory can be formulated as a bigravity mod- 
ification of General Relativity in which the gravity sector is extended by a new 
metric- like 2-tensor that interacts with the original metric only in particular ways. 
They also showed that as long as the structure of the theory is preserved, the new 
field can either be kept non-dynamical or be promoted to a dynamical one. 

Not only did the discovery of the ghost-free massive gravity, and its formula- 
tion in terms of a bimetric theory, attract the attention of many theoreticians, it 
also introduced a consistent, very interesting and theoretically well-motivated mod- 
ification of gravity to the cosmological community (for a review of cosmological 
applications of modified gravity theories, see Ref. 11). Such a non-linear theory of 
massive gravity may in particular explain cosmic acceleration without introducing 
the problematic cosmological constant (CC)^^ or Dark Energy. ^'^ 

In a longer paper, we investigated the ability of the bigravity theory of massive 
gravity to address the observed acceleration of the Universe^^'^^ as an alternative 
to the ACDM standard model of cosmology. We summarize those findings here; the 
reader is referred to the longer paper for details of all calculations, numerical and 
statistical methods used and detailed discussion of the results. 

We compare the background predictions of a cosmological model based on the 
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Table 1 . Best-fit > p- value and log-evidence for the ghost-free bigravity model and its sub- models when 
constrained by various cosmological measurements at the background level. For each model, parameters 
that are allowed to vary are marked as "free" ; the non- varying parameters are fixed to zero. 
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Hassan-Rosen bigravity theory to different cosmological observations. The two met- 
rics of the model are assumed to be spatially flat, homogeneous and isotropic. Only 
one metric is coupled to matter and is considered as the physical metric that ap- 
pears in the definitions of all cosmologically interesting observables. The full version 
of the model possesses six free parameters: il'^^, the present value of the (dark plus 
baryonic) matter density parameter, and Bi (z = 0, ...,4), representing possible in- 
teractions between the two metrics (see Refs. 14 and 17 for details). Bq corresponds 
to the explicit CC for the physical metric. We are therefore mostly interested in sub- 
models where Bq is set to zero (in order to test the ability of the model to provide a 
self-acceleration mechanism) . We perform an extensive Bayesian parameter estima- 
tion procedure using MultiNest algorithm^® and compute uncertainties around the 
best-fit points using marginalized posterior probabilities. We in addition employ a 
semi-frequentist technique to verify whether the parameters are well-constrained or 
there are correlations between them. We also compare the model's ability to fit the 
data to that of ACDM by evaluating Xmm ' P- value and log-evidence for each model. 
The data we use include the position of the first peak on the cosmic microwave back- 
ground angular power spectrum, the ratio of the sound horizon at the drag epoch 
to the dilation scale at six different redshifts,^°"^^ the luminosity distances to 580 
Type la Supernovae,^"^ and the present value of the Hubble parameter. 

Our numerical results, confirmed with detailed analytical studies, show (see Ta- 
bled]) that the model in general gives very good fits to the data at a very high confi- 
dence level. Except for (B2, ^Ij^), (B3, ^°m), {B2, B4, O^) and {B3, B4, n^,) that are 
ruled out, other sub-models with Bq — are in perfect agreement with observations 
and are statistically as consistent with the data as ACDM (log-evidence should 
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not be considered here because of the prior choices of parameter ranges and prior 
preferences we may have for a model over the other). In addition, we observe that 
models with more than one B are degenerate, i.e. parameters are correlated. Such 

degeneracies can be broken by using more data in a perturbativc analysis. Only for 
f2^), where we have reliable constraints, the graviton mass can be determined; 
it is of the order of the present value of the Hubble parameter. 

Acknowledgements 

We thank Tessa Baker, Robert Crittenden, Jonas Enander, Hans Kristian K. Erik- 
sen, Farhan Feroz, Pedro G. Fcrreira, Juan Garcia-Bellido, S. F. Hassan, Antony 
Lewis, Edvard Mortsell, David F. Mota, Sigurd K. Naess, Claudia de Rham, Rachel 
A. Rosen, Mikael von Strauss and Andrew J. ToUey for enlightening and helpful 
discussions. YA is supported by the European Research Council (ERC) Starting 
Grant StG2010-257080. 



References 

1. C. de Rham and G. Gabadadze, Phys.Rev. D82, p. 044020 (2010). 

2. C. de Rham, G. Gabadadze and A. J. Tolley, Phys. Rev. Lett. 106, p. 231101 (2011). 

3. S. F. Hassan and R. A. Rosen, JHEP 07, p. 009 (2011). 

4. S. Hassan and R. A. Rosen, Phys. Rev. Lett. 108, p. 041101 (2012). 

5. S. F. Hassan and R. A. Rosen, JHEP 02, p. 126 (2012). 

6. S. F. Hassan, R. A. Rosen and A. Schmidt-May, JHEP 02, p. 026 (2012). 

7. S. F. Hassan and R. A. Rosen. JHEP 04, p. 123 (2012). 

8. M. Fierz and W. Pauli, Proceedings of the Royal Society of London. Series A, Math- 
ematical and Physical Sciences , 211 (1939). 

9. S. Deser and D. G. Boulware, Phys. Rev. D6, 3368 (1972). 

10. K. Hinterbichler, Rev. Mod. Phys. 84, 671 (2012). 

11. T. Clifton, R G. Ferreira, A. Padilla and C. Skordis, Phys.Rept. 513, 1 (2012). 

12. J. Martin, Comptes Rendus Physique 13, 566 (2012). 

13. E. J. Copeland, M. Sami and S. Tsujikawa, Int.J.Mod.Phys. D15, 1753 (2006). 

14. Y. Akrami, T. S. Koivisto and M. Sandstad (2012). 

15. A. G. Riess et al, Astron.J. 116, 1009 (1998). 

16. S. Perlmutter et al., Astrophys.J. 517, 565 (1999). 

17. M. V. Strauss, A. Schmidt-May, J. Enander, E. Mortsell and S. F. Hassan, JCAP 03, 
p. 042 (2012). 

18. F. Feroz, M. Hobson and M. Bridges, Mon.Not.Roy.Astron.Soc. 398, 1601 (2009). 

19. E. Komatsu et al., Astrophys.J. Suppl. 192, p. 18 (2011). 

20. F. Rentier, C. Blake, M. CoUcss, D. H. Jones, L. Staveley-Smith et al, 
Mon.Not.Roy.Astron.Soc. 416, 3017 (2011). 

21. W. J. Percival et al, Mon.Not.Roy.Astron.Soc. 401, 2148 (2010). 

22. C. Blake, E. Kazin, F. Beutler, T. Davis, D. Parkinson et al., Mon.Not.Roy.Astron.Soc. 
418, 1707 (2011). 

23. N. Suzuki, D. Rubin, C. Lidman, G. Aldering, R. AmanuUah et al., Astrophys.J. 746, 
p. 85 (2012). 



